Bound states of PT-symmetric separable potentials 
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All of the PT-symmetric potentials that have been studied so far have been local. 

In this paper nonlocal 'PT-symmetric separable potentials of the form V{x, y) = 
^~~' ie\U{x)U{y) — U{—x)U(—y)], where U{x) is real, are examined. Two specific models 

(^ are examined. In each case it is shown that there is a parametric region of the 

CN coupling strength e for which the VT symmetry of the Hamiltonian is unbroken and 

the bound-state energies are real. The critical values of e that bound this region are 

calculated. 
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I. INTRODUCTION 

The time-independent Schrodinger equation Hip = Eip, when written in coordinate space, 
takes the form of a continuous matrix eigenvalue problem 



Oh 

> 

m r 

0> dyHix,y)iPiy) = E,Pix). (1) 

(N 

(^ However, all previous studies of non-Hermitian PT-symmetric quantum-mechanical Hamil- 

O tonians have focused on Hamiltonians that in the coordinate representation are diagonal 

I and symmetric. Such Hamiltonians take the form H = p"^ + V{x), where the condition of 

IL" VT symmetry is that V*{x) = V{—x). When expressed as a matrix, H is clearly diagonal, 

>< Hix, y) = d,dy5{x -y) + V{x)5{x - y), (2) 

and for Hamiltonians of this form the Schrodinger eigenvalue problem (fTl) is the differential 
equation 

- ^"{x) + v{x)ij{x) = E^{x). (3) 

The first PT-symmetric Hamiltonians that were examined in detail belong to the class 
of diagonal Hamiltonians P-|l] 

H = p^ + x^{ixy (e>0). (4) 

The Hamiltonians Q can have many different spectra depending on the large- |x| bound- 
ary conditions that are imposed on the solutions to the corresponding time-independent 
Schrodinger eigenvalue equation 

- i)'\x) + x'^{ixyi){x) = Eip{x). (5) 
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The boundary conditions on iIj{x) are imposed in Stokes' wedges in the complex-x plane. 
At the edges of the Stokes' wedges both hnearly independent solutions to (tsl) are oscillatory 
as |a:| — >■ oo. However, in the interior of the wedges one solution decays exponentially and 
the linearly independent solution grows exponentially. The eigenvalues E are determined 
by requiring that ip{x) decay exponentially in two nonadjacent wedges. Ordinarily, the 
eigenvalues are complex, but if the two wedges are "PT-symmetric relative to the imaginary- 
X axis, then all the eigenvalues are real. (The VT reflection of the complex number x is the 
number —x*.) 

We emphasize that a Hamiltonian need not be diagonal. Of course, for nondiagonal 
Hamiltonians it is difficult to solve the Schrodinger eigenvalue problem because it takes the 
form of the Volterra intego-differential equation ([I]) rather than the differential equation (Is]). 
However, there is an interesting special solvable class of nondiagonal Hamiltonians for which 
the potential is separable; these Hamiltonians have the form 

H{x, y) = d,dy5{x -y) + U{x)U{y). (6) 

Potentials of this form are discussed in Ref. [5j and are interesting because they can be used 
in studies of scattering processes. 

The purpose of this paper is to consider the case of complex PT-symmetric separable 
potentials. We show that such potentials can have unbroken (and broken) PT-symmetric 
regions in which the bound states are real (and complex). In Sec. Hlwe show how to find 



bound states for real Hermitian separable potentials and in Sees. HI and IV we show how to 
find bound states for complex PT-symmetric potentials. Finally, in Sec. |V] we make some 
brief concluding remarks. 



II. BOUND STATES OF HERMITIAN SEPARABLE POTENTIALS 

Let us consider the case of a Hamiltonian of the form H = p^ + V{x,y), where V{x,y) = 
gU{x)U{y) is a separable potential and (yf is a coupling strength. Note that if U{x) is real, 
then H is Hermitian. The Schrodinger eigenvalue equation for this Hamiltonian is 

- r{x) + gU{x) / dyi^{y)U{y) = E^{x), (7) 



which is a linear integro-differential equation in Volterra form. 

In general, the potential V{x,y), may have bound states and scattering states. These 
states are distinguished by their large-|x| behavior. For a scattering state the wave function 
ip{x) does not vanish as \x\ -^ cx), but for a bound state ipi^) — )■ as \x\ — >■ oo. Furthermore, 
if ip{x) vanishes for large |a;|, then ([T]) implies that U{x) must also vanish. Thus, if we wish 
to solve ([T]) for a bound state, it is valid to perform a Fourier transform: 



{E-p')^{p)=gU{p) dy^{y)U{y), (8) 

J — oo 

where the Fourier transform of %l){x) is given by 

i){p)= I t/xe^^Xx) (9) 



and the inverse Fourier transform of ipip) is given by 

^{x) = — / dpe-'P''^ip). (10) 

2vr J-oc 

Note also that the integral in (|8| can be written as a convolution of two Fourier transforms: 

dy^{y)U{y) = — / dqU{-q)i,{q). (11) 

To proceed, we multiply ([8| by U{—p) = U*ip), integrate with respect to p, and assume 
that 

/oo 
dqUi-q)i,{q)^0. (12) 

-oo 

The result is the secular eigenvalue condition 



2n 



dp 



\uip)? 

E-p"^ 



(13) 



The integral in this equation must exist, and it clearly exists if E has a nonzero imaginary 
part. However, complex energy is associated with scattering, and we are not concerned with 
scattering states in this paper. 



There are two possible ways to have real-energy bound-state solutions to (13). First, E 



may be negative. In this case the integral exists but it is negative, and this requires that the 
coupling constant g be negative. While this possibility is viable for the case of Hermitian 
Hamiltonians, it is not viable for the case of non-Hermitian PT-symmetric separable Hamil- 
tonians, as we will show in Sec. ul Another way to have a bound state is for U{p) to vanish 
when p^ = E, where E is a number to be determined consistently, and if this happens, the 



integral in (13) will exist and there can be a bound state. 



Let us construct a potential for which we can solve the one-bound-state problem analyt- 
ically: We assume for simplicity that U{p) has the form 



U{p) 



-p2/2 



{E-p' 



Then, from evaluating the integral in (13) we get an equation for E in terms of g: 



We also obtain the result that 

U{x) = 

We now verify the consistency o 
To do so, we calculate ip{p) from 



E 



1 

/2^ 



1 20^ 
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-2+"' 



-x'^/2 



(14) 



(15) 



(16) 



ga{2TT) 



-3/2 -xV2 



this calculation by showing that a in (12) is nonzero 



8j) and find that ip{p) = gae ^ '^/(27r), and thus that 



We then evaluate (12) and find that it reduces to the identity 



a = a. Thus, a is an arbitrary normalization constant and the bound-state solution is 
internally consistent. 



III. BOUND STATES OF PT-SYMMETRIC SEPARABLE POTENTIALS 



For the non-Hermitian case we assume that the potential has the PT-symmetric form 
V{x,y) = ei[W{x,y) — W{—x,—y)], where W{x,y) is reaL We then further speciahze 
the potential by making it separable and symmetric under the interchange of x and y: 
V{x,y) = ie[U{x)U{y) — U{—x)U{—y)], where U{x) is real. This gives the Schrodinger 
equation 



ijj"{x) +ie 



We can rewrite this equation as 



f/(x) / dyi;iy)Uiy)-Ui-x) / dyi^iy^i-y) 



le 



'^'\x) + 7^ W{x) - l3U{-x)] = Eij{x), 



Ei){x). (17) 



(18) 



where a and /3 are expressed as convolutions of the Fourier transform of U and the Fourier 
transform of ip: 



a 



dqdhqMq), /3 



dqU{q)i}{q). 



(19) 



We assume that a and /3 are nonzero. 

As in the Hermitian case discussed in Sec. [Ill we solve the Schrodinger equation (18) by 
taking a Fourier transform: 



^{p) 



2^ 



u{p) uj-p) , 

E — p^ L — p^ 



(20) 



Then, we multiply (20) by U{—p) and integrate to obtain 

mp)\' 



a 



2^ 



a 



dp 



E-p^ 



/3 



dp 



and we multiply (20) by U{p) and integrate to obtain 

[U{p)f 



P 



le 

2^ 



a 



dp 



E-p^ 



(3 / dp 



E-p"^ 



E-p'^ 



(21) 



(22) 



These two equations have the form 



where 



a 



h 



S"'^ - ^^'^' 



/3 






dp 



{u{p)f 

E-p^ 



dp 



\U{P)? 
E-p^ 



(23) 



(24) 



Note that U{p) must be complex [otherwise (23) reduces to triviality a = /3 = 0]. It is 
necessary to assume here that /i and I2 exist. This requires that U{p) vanish at p^ = E. 
Thus, the eigenvalues of the Hamiltonian must be zeros of U{p). (Other solutions for which 



E is complex correspond to scattering states. As stated earlier, in this paper we consider 
only bound states.) 



Because the simultaneous linear equations (23) are homogeneous a solution for a and f3 



exists only if the determinant of the coefficients vanishes. This requirement gives a secular 
equation that determines the energy E: 



e 

47^2 



i^i - ^2) 



(25) 



which is the analog of ( 13 ) for the case of a Hermitian separable potential. 



This secular equation has several important properties and consequences. First, like the 
secular equation for all "PT-symmetric Hamiltonians, it is real if E is real [B]. Note that I2 in 
(24) is manifestly real. To see that /i in ( |24[ ) is real, observe that because U{x) is real, U{p) 
is VT symmetric; that is, U*{p) = U{—p). Thus, the change of variable p — )■ —p establishes 
the reality of Ji. As a consequence, the roots E of this equation are either real or come in 
complex-conjugate pairs. Second, (25) in conjunction with either of the two equations in 



(23) implies that 



l«l = l/3|- (26) 

Third, observe that £" is a function of the square of the coupling constant e so it is inde- 
pendent of the sign of e. This feature is typical of PT-symmetric Hamiltonians such as 
H = p'^ + x"^ + iex and H = p^ + x"^ + eix^, which have imaginary potentials. Fourth, there 



is no solution to (25) if E is negative. This is because when E < the triangle inequality 
implies that /| > I{, and thus the sign of the right side of (25) is negative for any choice of 
e. Thus, if there is a bound state, its energy must be positive. _ 

To have a bound state for positive energy E, the function U{p) must vanish at p"^ = E. 
As we did in the Hermitian case discussed in Sec. [TT| we construct a simple model for which 
there is just one real root of U{p) at E = p^: 



U{p) 



V/2 



(E -/)(! + zap). 



where a is a real constant so that U{x) is a real function, namely, 

1 



U{x) 



^-x ,1^^ _ 1 + (^ _ 3)^^ + x' + ax' 



(27) 



(28) 



[Note that if we had chosen this U{x) first, it would not have been obvious that U{p) had a 
zero.] Then 

I^ = ^(e--- -a^E + -a^ 
^ ^ \ 2 2 4 



TT \ E h -a^E a^ 

2 2 4 



Hence, 



h-h 



1 



'V^{3-2E), h + l2 = V^{2E~l). 



and the secular equation (25) reduces to the quadratic equation 

+ {2E - 1){2E - 3) = 0, 



a^e^ 



(29) 



(30) 



(31) 
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FIG. 1: The bound-state energy Em. (32) plotted as a function of ae. The region of unbroken VT 
symmetry is ae > yovr = 5.013 .... In this region the Hamiltonian has only one bound state, but 



the two allowed values of the energy approach the asymptotic limits I and ^ as ae 



there are two possible allowed values for the energy of this bound-state. As one can see from (32), 

1 

2 "^^^ 2 
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whose roots are 



E = 1± 



1- 



iir 



a^t^ 



(32) 



This equation gives a condition on the coupling constant e for the reality of the bound-state 
energy E: 

e > ^. (33) 

a 

The interpretation of this inequality is straightforward. Even though U{x) — )■ for large 

|a;| it is still possible for there to be a bound state if e is large enough. However, when e 

is smaller than a critical value, the bound state disappears (its energy becomes complex). 

This corresponds to going from a PT-unbroken region to a PT-broken region. A plot of E 

as a function of ae is given in Fig. [ij Note that the Hamiltonian has only one bound state; 

however, there are two allowed values for the energy of this bound state. 

Finally, we must verify that apart from a constant multiplicative phase [which is arbitrary 

because the Schrodinger equation (17) is linear] the eigenfunction ip{x) is VT symmetric. 

To do so, we calculate ip{p) and verify that apart from a constant multiplicative phase the 

function il){p) is a real function of p. To obtain ^(p) we substitute (27) into (20) and obtain 



V'(p) 



^ e-p''^[i{a-(5)-ap{a + (3)] 



271 

e 



-(«.«e-.'/^(,^ 



ap 



(34) 



Next, we add the two equations in (23) and immediately obtain the result that the ratio 
[a — (3)/{a + (3) = —27ii/{Ii + I2) is imaginary. This verifies that up to an arbitrary 
multiplicative phase, iplx) is VT symmetric; it becomes VT symmetric when a + /3 is 
chosen to be real, so that a = (3*. If the coupling constant e lies below the critical point 



7 



e < 27r-\/2/a, then E becomes complex; consequently, /i and I2 are not real and tp{x) is not 
VT symmetric. 



IV. SEPARATED PT- SYMMETRIC POTENTIAL HAVING TWO BOUND 

STATES 

Although it becomes technically complicated, it is straightforward to generalize the dis- 



cussion of Sec. |III| to the case in which the PT-symmetric separated potential has more than 
one bound state. In this section we consider the case in which there are two bound states 



of energy Ei and E2 for the Schrodinger equation (17). 



The bound-state eigenfunctions are labeled ipkix) {k = 1,2), and in momentum space 

U{-p) 



they have the same general form as in ( 20 ) : 

U{p) 



^k{p) 



le 

2^ 



Eh 



pz 



-,otk 



Eh 



pZ 



;Pk 



(A; = 1,2), 



where 



ak 



dqU{-q)iJk{q), (3k 



dqU{q)Mq) (fc = l,2). 



(35) 



(36) 



Equations (|21|) and (|22|) generalize to 

|/7(p)p 



Otk 



h 



le 

2^ 

ie 

2^ 



otk 



Oik 



dp 
dp 



Ek-P' 

[UJP)? 
Ek-p^ 



h 



(3k 



dp 
dp 



[u{p)f 

Ek~p' 

mp)\' 

Ek-p^ 



(A; = 1,2), 
(fc = l,2). 



(37) 



We then extend (27) so that U{p) has two real zeros instead of one: 

U{p) = e-'P"/\E,-p^){E2-p^){l + lap). (38) 

The secular equations satisfied by the energies Ei and E2 are a rather complicated general- 



ization of (31): 



1287r 

a?e^ 
1287r 



+ P{E,,E2)Q{E,,E2) = 0, 
+ P{E2,Ei)Q{E2,Ei) = 0, 



(39) 



where 



P{E,,E2) 
Q{E,,E2) 



-15 + 6E1 + I2E2 - 8E1E2 - 4:El + 8EiEl, 
-105 + 30^1 + 6OE2 - 24^1^2 - l2El + SE^E. 



2- 



(40) 



By solving simultaneously the equations (39) and (40), we obtain the allowed bound state 



energies. The numerical solution of these equations is plotted in Fig. [2j There is a critical 
point at ae = 1.09; below this value there are no real solutions. Bound states first appear 
when ae exceeds this critical value and there are two cases: (i) There can be just one bound 
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FIG. 2: Numerical solution to the simultaneous equations (39) and (40). Bound-state energies 



are graphed as functions of ae. The region of unbroken VT symmetry is ae > 1.09, and bound 
states appear as soon as ae exceeds this critical value. A second critical point is at ae = 3.90, at 
which four new solutions appear. There are two cases: In the first case there is only one bound 
state whose energy may lie on the curve d or on the curve e. In the second case there are two 
bound-state energies, which may lie on the curves a and /, or on b and g, or on c and h. The two 
curves g and h are too close to be resolved in this figure and are thus shown in detail in Fig. [3} As 
ae — )• oo, the curves a - h approach the asymptotic values 4.081, 3.742, 2.725, 2.115, 1.054, 0.919, 
0.296, and 0.275. 

state if El = E2, and in this case there are two possible values for this energy and these are 
indicated in Fig. |2]by the curves consisting of connected dots and labeled d and e. (ii) There 
can be two distinct bound states; these are indicated by the sequences of dots labeled a and 
/. If El lies on the curve a, then E2 lies on /, and vice versa. It is interesting that there is a 
second critical point at ae = 3.90 and when ae exceeds this value four new possible energies 
appear; these energies lie on the curves labeled by b, c, g, and h. The possible bound-state 
energies Ei and E2 are only allowed to lie on the pairs of curves a and /, or b and g, or c 
and h. As ae — )■ 00, the energies on the curves a - h approach the asymptotic values 4.081, 
3.742, 2.725, 2.115, 1.054, 0.919, 0.296, and 0.275. The energies on the curves g and h are 
too close together to be resolved in Fig. |2| so we have included a separate figure Fig. |3] to 
show their dependence on ae. 



V. COMMENTS AND DISCUSSION 

We have shown in this paper that it is easy to construct non-Hermitian PT-symmetric 
separable potentials and, even though such potentials are nonlocal, it is still possible to find 
the secular equation that determines the bound-state energies. As is the case for any VT- 
symmetric potential, the secular equation is real. If one solves the secular equation, one finds 
a result that is typical of PT-symmetric theories; namely, that the coupling constant lies in 
one of two regions, which are separated by a critical value: On one side of the critical value 
(the unbroken region of unbroken VT symmetry) the energies are real, but on the other side 
of the critical point (the region of broken VT symmetry) the energies are complex. This VT 
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FIG. 3: Blow-up of the curves g and h from Fig. [2] 

phase transition has been observed experimentally in PT-symmetric optical models [TIIH]. 

Based on the structure and behavior of the models we have constructed in this paper, it 
is evident that we can construct PT-symmetric separable potentials for which there are as 
many bound states as we wish, and in the unbroken PT-symmetric region all of the bound 
states will have positive energies. 

For the models discussed in this paper the potentials vanish exponentially rapidly as 
|x| — )■ oo. Thus, for large |x| the Hamiltonian becomes the free Hamiltonian Hq = p"^, 
whose solutions are plane waves. Thus, in addition to bound states, there will be scattering 
states for all positive energies. The energy of a scattering state will be complex, with the 
sign of the imaginary part of the energy being associated with incoming- or outgoing-wave 
boundary conditions. The models we have studied in this paper are interesting because the 
point spectrum of bound states is embedded in the continuum of scattering states [9HTT]. 
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